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We explore the non-equilibrium dynamics of non-interacting classical particles in a one- 
dimensional driven superlattice which is composed of domains exposed to different time-dependent 
forces. It is shown how the combination of directed transport and conversion processes from diffu- 
sive to ballistic motion causes strong correlations between velocity and phase for particles passing 
through a superlattice. A detailed understanding of the underlying mechanism allows us to tune 
the resulting velocity distributions at distinguished points in the superlattice by means of local vari- 
ations of the applied driving force. As an intriguing application we present a scheme how initially 
diffusive particles can be transformed into a monoenergetic pulsed particle beam whose parameters 
such as its energy can be varied. 
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I. INTRODUCTION 

Dynamical systems and their transport properties 
have been studied extensively over the last decades [l|, 
whereas the topic stimulated a vast amount of research 
when the possibility of directed currents in the absence 
of a mean force was realised 043 ■ Since the second law 
of thermodynamics forbids such transport phenomena in 
equilibrium, these systems have to be constantly driven 
out of equilibrium. Early works [2hJ] were based on noise, 
i.e. statistical external fields, in combination with spa- 
tially asymmetric so called 'ratchet' potentials to over- 
come the limitations formulated by the second law of 
thermodynamics and thus to evoke a particle current. 
These type of systems are of particular interest because 
they outline a working principal for biological systems 
such as molecular motors [5|, or quantum motors pj . 
However, it was soon realised that directed currents can 
very well be obtained with deterministic external fields, 
as long as certain spatial- and temporal symmetries in 
the equations of motion are broken |8J, which was in- 
vestigated afterwards in a vast amount of literature (see 
p-fla] and references therein). These deterministic ratch- 
ets are of particular interest since they might have re- 
markable applications in nanoscale devices such as elec- 
tron pumps or transistors 16|. First experimental reali- 
sations included semiconductors or semiconductor micro 
structures where a combination of Laser fields has been 
applied which led to directed currents in electron ratch- 
ets [13, [l8|. Directed currents also became a subject of 
interest in experiments concerning cold atoms in optical 
lattices [19l - l22| , where additional AC forces are applied to 
drive the system out of equilibrium and at the same time 
break the required symmetries that would otherwise pre- 
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vent transport phenomena. These type of experiments 
are of particular interest since they allow for a precise 
control over the system parameters and provide exten- 
sive tuneability in the used AC drivings [20, [22[ . 
While the so far mentioned works focus on only time- 
dependent AC forces, it was found recently that a spa- 
tial dependence of these AC forces leads to a diversity 
of dynamical phenomena 23-25]. The latter studies ad- 
dress the classical dynamics of particles in a lattice with 
a site-dependent driving. In [23| it is demonstrated how 
a phase-modulated lattice allows for directed transport 
even though the driving of each barrier on its own does 
not break the relevant symmetries [8(. Ref. [2J] shows 
how a ramping of the potential height in combination 
with a site-dependent driving leads to a patterned de- 
position of particles. A specific realisation of a site- 
dependently driven lattice is the one of a block lattice 
(BL) as introduced in [25|], which is reminiscent of semi- 
conductor heterostructures and superlattices. Indeed 
only recently the possibility of ratchet effects in super- 
lattices of semiconductor heterostructures with a super- 
imposed periodic potential was reported [26|, [27| . In the 
case of [25], the superlattice consists of different blocks 
containing many inidvidual barriers where the barriers of 
each block are governed by a certain time-dependence i.e. 
driving law, whereas different blocks exhibit in general 
different driving laws. The long time transient dynamics 
in such a superlattice shows intriguing phenomena like 
the formation of spatial density oscillations. The latter 
were explained and analysed by means of conversion pro- 
cesses fronr diffusive to ballistic (and vice versa) motion 
at the positions where two neighbouring blocks connect, 
i.e. at the interfaces of two blocks. However, a rigorous 
discussion of the processes occuring in a unit cell of such a 
superlattice -that is a system containing only two blocks 
each equipped with one of the used driving laws- is still 
missing and is therefore subject of the present work. In 
this sense we investigate the diffusive- to ballistic motion 
conversion processes in detail and explore their influence 



on the dynamics of particles leaving the two block sys- 
tem. As a result we obtain peaked velocity distributions 
for outgoing particles even though their initial conditions 
are chosen exclusively within the chaotic sea of the under- 
lying phase space. By adjusting parameters in the driving 
we are able to manipulate these velocity distributions in 
a controlled manner. Finally, we demonstrate how the 
insights gained from the two block system enable us to 
exploit the conversion processes in superlattices build up 
of many blocks each equipped with a unique driving law. 
In doing so an initially diffusive particle ensembles is con- 
verted into an ensemble with a velocity distribution con- 
taining a single peak, whereas both the width as well as 
the peaks mean velocity can be tuned. 
The present work is structured in the following way: In 
section QI] we introduce the setup of a block lattice (BL). 
In section Uni the dynamics of a single block is discussed. 
We explore the conversion of diffusive to ballistic mo- 
tion at distinguished positions in the BL in section IIVI 
Additionally, the influence of these processes for outgo- 
ing particles in a simple two block system is discussed in 
section [Vj Finally, we investigate the dynamics of super- 
lattices containing several hundred blocks in section |VT1 
Section [VIII contains our brief conclusions. 



II. SETUP AND HAMILTONIAN 

The system investigated is a one-dimensional driven 
lattice consisting of laterally oscillating square potential 
barriers of equal height Vq and length I as sketched in Fig. 
[TJ Each barrier is characterised uniquely by its index i. 
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FIG. 1. Sketch of a lattice that consists of differently laterally 
driven blocks, where the barriers with NB < i < (B + 1)N 
are equipped with the same driving law ds(i). 

Furthermore, the barrier positions are time-dependent 
and described by the so-called 'driving law' d(t), which 
is chosen such that the i-th barrier remains at all times 
within an interval of length L expanding from iL to 
(i + 1)L. Additionally, the lattice is divided into blocks 
such that different driving laws dsit) (introducing the 
'block index' B) are used. In doing so, each dsit) gov- 
erns the barrier motions for the sites NB < i < (B + 1)N 
where N denotes the number of barriers within one block 
and is set to 10 4 throughout this work. The general 
structure of the driving law is a biharmonic function 
rfs(t) = AB[cos(ujBt) + sin(2ujBt + A$b)] with three pa- 
rameters Ab, i^b an d A$b which depend on the block 
index B. Hence, the resulting classical Hamiltonian for 
noninteracting particles is given by 



B=N B1 -1(B+1)N (1) 

J2 E v e(!/2-\x-Xo,i-d B (t)\), 

B=0 i=BN 

with Abi being the number of considered 
blocks and Xq^ the equilibrium position of the 
j-th barrier (chosen such that the barrier os- 
cillates symmetrically within its unit cell, i.e. 
min|X 0j j + d B {t) - iL\ = min|A 0:l + d B {t) - (i + 1)L\). 
Additionally, we set the mass m — 1 without loss of 
generality and keep Vb = 1.0, L — 5.0 and I = 1.0 
constant throughout this work. 



III. DYNAMICS IN A SINGLE BLOCK 

Even though the focus of this work is on compos- 
ite systems consisting of multiple blocks exposed to 
different driving laws, these blocks are considered to 
be large in a sense that the dynamics of a particle 
within one block can be described by the Poincare 
surfaces of section (PSS) as obtained by extending 
this block to an infinite uniformly driven lattice. It 
is therefore sensible to discuss the dynamical features 
such as the transport properties as well as the appear- 
ance of the PSS of the uniformly driven lattice. It is 
well established that the PSS for a uniformly driven 
lattice can be obtained by exploiting the temporal 
symmetry of the Hamiltonian: H(x,t) — H(x,t + T). 
Thus an appropriate choice for the PSS is given by 
M = {{x(t + kT) mod L, p{t + kT)) \ k G N} with 
T = ^- being the temporal period. According to this, 
we record the position taken modulo L and the velocity 
at certain times 2irk, k G N (and call this 'position 
velocity section'). Such a PSS for cj = 1.0, A « 0.57 and 
A$ = is shown in Fig. [2] b) and features the typical 
mixed phase space [28(, i.e. KAM islands embedded in 
a chaotic sea which is bounded by the first invariant 
spanning curve (FISC). Note that the white rectangle is 
caused by adding the potential energy for particles which 
are within a barrier at times when position and velocity 
are recorded. This is done to avoid discontinuities in the 
PSS caused by the discontinuous potential V(x, t) (cf. 
[231]). For later usage we show additionally the position 
velocity section for parameters: uj = 1.0, A w —0.57 and 
A$ = in Fig. Ha). Note that Fig. Ha) (b)) shows the 
PSS of the left block (right block) of the corresponding 
two block system (see section HV| . 

A second possibility to illustrate the phase space of a 
uniformly driven lattice is given by the 'phase velocity 
section' which exploits the spatial symmetry of the 
Hamiltonian: H(x,t) — H(x + L,t). To this end we 
record the phases and velocities at certain positions 
iL. Hence the PSS is obtained from the set of points 



M = {{t(x + kL) mod T, p(x + kL)) \ k € N} and the 
resulting plot is shown in Fig. [3] d) (again, the PSS for 
parameters as in Fig. [2]a) is shown additionally in Fig. [2] 
c) for later usage). Apparently, it features qualitatively 
the same domains as the position velocity section (Fig 
[2] b)), i.e. ballistic islands which are embedded in a 
bounded chaotic sea. However, in contrast to the pre- 
vious case the chaotic sea appears to be non uniformly 
filled with trajectories. This seeming contradiction to 
ergodicity can be resolved easily: According to ergodicity 
the chaotic sea in Fig. f2]b) can assumed to be filled with 
a uniform measure. Hence the number of particles AN 
that pass Xi ~ iL and therefore contribute to the phase 
velocity section per time At and velocity interval Av 
is given by -^j- = p^jAv — pvAv, where Ax denotes 
the distance a particle travels in time At and p is the 
number of particles per phase space interval. Therefore, 
the number of particles passing Xi per velocity interval 
is ^- — pvAt and hence proportional to v. 
We now comment on the transport properties within 



Fig. [3] and appear to be in good agreement to the 
results of a symmetry analysis (cf. [g, [29|) which yields 
^transport oc — cos(A<f>). However, there are noticeable 
deviations, e.g. the reversed sign close to A$s = n/2 
and A$s = 3ir/2. Note that these deviations should 
not surprise us because the authors in [3, |29| considered 
continuous potentials instead of discontinuous potential 
barriers. For a more detailed analysis of the dynamics 
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FIG. 3. Numerically determined transport velocity as a 
function of A$ in an infinite uniformly driven lattice with 
A w 0.57 and u = 1.0. 



in a uniformly driven lattice, we refer to [23j | where this 
was done in great detail. 



IV. INTERFACE CONVERSION IN THE TWO 
BLOCK SYSTEM 

In Ref. |25| it was argued that a BL as introduced in 
section |TT] offers the opportunity for conversion processes 
from diffusive- to ballistic motion and vice versa at 
interface positions where the driving law changes. These 
type of processes will be analysed in detail throughout 
this section. We demonstrate in particular their in- 
fluence on the velocity distribution of a particle ensemble. 



FIG. 2. a) and c) PSS for a uniformly driven lattice with 
parameters A « 0.57, A$ = and u> - 1.0. b) and d) PSS 
for A m —0.57, A$ = and u = 1.0. a) and b) phases and 
velocities are recorded at positions x mod L = 0. c) and d) 
position x taken modulo L as well as the velocity at times 
tot mod 2ty = 0. 

a single block. For A$s ^ n\ (n E Z) the biharmonic 
driving law breaks the time-reversal invariance as well as 
the parity symmetry of the Hamiltonian and the driven 
lattice allows for directed transport phenomena (8]. The 
transport as a function of A<& with fixed A « 0.57 and 
lu = 1.0 is determined numerically by simulating 10 
particles in a uniformly driven lattice for 10 6 barrier 
oscillations and calculating their average velocity after 
a certain transient time. The results are shown in 



A. Interface conversion 

Let us introduce the simplest possible finite BL which 
is build up of only two blocks (i.e. JVbi = 2) equipped 
with different driving laws. Such a system extends from 
^min = to X max = 2N L (so the simulation is stopped for 
a particle once it passes either of these positions) and the 
driving laws are do(t) for x < x nll d — NL and d\(t) for 
x > iTmid- The parameters of di(t) are chosen as before 
(wi = 1.0, A\ « 0.57 and A$i = 0) and thus the dynam- 
ics within the 'right block' (RB) can be described by the 
two PSS in Fig. [5]b) and d). Moreover, Fig. [3] reveals 
that the used driving law induces a negatively directed 
current. For the 'left block' (LB) we chose ojq = 1.0, 



A sa -0.57 and A$ = yielding d Q (t) = -di(t). Hence 
the corresponding position-velocity section is given by 
Fig. [2] a) and the phase- velocity section is the one shown 
in Fig. d]c). The induced current in the LB is therefore 
of the same magnitude as in the RB, but positively di- 
rected. 

To understand how this setup allows for conversion pro- 
cesses it is helpful to consider the dynamics of a parti- 
cle with initial conditions in the chaotic sea of the LBs 
phase space. Due to the positively directed current, this 
particle is in the average transported towards x m ^ and 
the chaotic sea for positive velocities in the correspond- 
ing phase velocity section (Fig. [2]c)) marks all possible 
phase space coordinates (vd,4>d) at which the particle 
can reach x m id diffusively, while the coordinates belong- 
ing to ballistic motion (i.e. within ballistic islands or 
regular spanning curves above the FISC) {vb,4>b) are 
prohibited. However, once the particle passes the inter- 
face at x m id its dynamics is no longer governed by the 
LBs phase space, but by the phase space of the RB, 
which is appropriately described by the PSS in Fig. [5] 
d). The crucial observation is that some of the coordi- 
nates (vdt^d) belonging to diffusive motion in the LB 
correspond to ballistic motion in the RB. This is best 
seen by means of a concrete example: Imagine the par- 
ticle reaches x mi d with (v = 0.8,0 = 37r/2), which is 
inside the chaotic sea of the LBs PSS (Fig. [5]c)). For 
x > £ m id the particles dynamics is described by the RBs 
PSS (Fig. [5]d)) where these coordinates correspond to 
a ballistic island. Hence, this initially diffusive particle 
would have become ballistic at the interface and we refer 
to this process as diffusive to ballistic motion conversion. 
Besides being injected into ballistic islands, the particles 
can equally well be injected into regular curves above the 
RBs FISC, because the FISC for positive velocities in the 
LB is at higher velocities as it is in the RB. To state a 
general rule, initially diffusive particles can be injected 
into every regular structure of the RBs PSS which has at 
least some 'overlap' with the chaotic sea of the LBs PSS. 
An example of a regular structure in which no injection 
can occur is the chain of ballistic islands at v ~ 3.2 in 
the RBs PSS (Fig. [Id)). These islands are 'covered' by 
a chain of larger islands at the same velocity in the LBs 
PSS (Fig. He). 



B. Density evolution in the two block setup 

After having discussed the process of interface conver- 
sion in BLs, we will explore their influence on the time 
evolution of the particle density in the following. 
To this end we propagate the dynamics of a particle en- 
semble in the two block system, which we introduced in 
section I IV Al As initial conditions we chose uniform dis- 
tributions for the particles positions as well as their veloc- 
ities with OANL < x ini < 0.6NL and -0.1 < w ini < 0.1 
respectively. Hence the particles are symmetrically dis- 
tributed around the LBs center and started in the chaotic 



sea of the phase space. Naively one might expect that 
due to the oppositely directed currents in the LB and the 
RB an accumulation of particles might happen at the in- 
terface at x m id. As we shall see in the following, this 
does not occur due to the previously introduced conver- 
sion processes. 

Fig. 0]shows snapshots of the normalised particle density 
at different times. Fig. 0] a) shows the particle density 
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FIG. 4. Particle density at different times in the two block 
setup (TVbi = 2) and TV = 10 4 yielding a; m id = 5-10 4 (indicated 
by the red dashed line) and :r max = 10°. The two driving 
laws are do(t) with ojo = 1.0, Aa « —0.57 and A"l>o = for 
x < x m id and d\(t) = — do(i) for x > a; m id. All densities are 
normalised to the initial particle number 10 5 . Snapshots of 
the particle density are shown for t = (a), 10 (b), 10 (c), 
1.8 • 10 5 (d), 3 • 10 5 (e) and 5.5 • 10 5 (f). 

p(x) for t — 0. For t > the particle distribution starts 
to spread and reaches a Gaussian like shape at t = 10 4 
(Fig. 0]b)). Afterwards the ensemble drifts in positive 
x-direction and once a sufficient amount of particles ar- 
rives at x m id = 5 • 10 4 (marked by the red dashed line in 
Fig. [J) a sharp decrease of the density emerges at this 
position (Figs. [4] c) and d)). This effect outlasts until 
t « 5.5- 10 5 (Figs. He) and f)). Finally at t « 10 6 , 
all particles have left the system at either x = x m ; n or 
x = i max and thus the density in the system is zero. 
The broadening of the peak within the LB can be ex- 
plained by underlying diffusion processes, because all 
particles are initially restricted to the chaotic sea. In 
fact, as long as the particles have not reached x = x m - l( i 
the ensemble is super-diffusive [23|], leading to a fast ex- 
pansion. On the one hand, the observed average drift of 



the ensemble in the LB is explained easily by the posi- 
tively directed transport in the LB. On the other hand, 
the fast density decrease is -as mentioned before- some- 
what counterintuitive. However, this effect is a straight- 
forward consequence of the in section IIV Al introduced 
conversion processes at aWd- According to our discus- 
sion, the initially diffusive particles can be injected from 
the chaotic sea of the LBs PSS into regular structures of 
the RBs PSS which leads to a fast ballistic motion away 
from x m id- If the particle remains diffusive, the directed 
transport brings it back to the interface and an injection 
from the chaotic sea of the RB to regular structures of 
the LB is possible. Hence, this process is repeated un- 
til an injection occurs and the particle leaves the system 
at i max (or Xmin) within a regular structure of the RB 
(or the LB). The fact that we do not observe a parti- 
cle accumulation at the interface demonstrates that the 
conversion process from diffusive to ballistic motion hap- 
pens on a sufficiently fast timescale to overcompensate 
accumulation effects caused by the directed currents. 
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C. Analysis of escaping particles in the two block 
setup 

In the present section we investigate the conversion 
processes from diffusive to ballistic motions further 
and illustrate an intriguing hallmark of these processes, 
namely their influence on the particles phase-velocity 
distribution. 

The starting point is the same two block setup as before 
and the initial conditions are chosen as before, too. 
However, instead of discussing the particle positions 
at certain times as we did in section IIVBI we now 
record the particle velocities and phases at distinguished 
positions: rr m i n and x max . To this end the phase (f> and 
the velocity v for every particle at x max are recorded 
and the result is shown in Fig. [5] a). In the low velocity 
regime (v < 3.0) a distinguished island structure is 
apparent, whereas for higher velocities (v > 3.5) the 
particles possess all possible phases from <fr — to 
4> = 2tt. In between (3.0 < v < 3.5) the particles appear 
to have randomly distributed phases and velocities, but 
do not occupy certain islands. 

The islands in the low velocity regime are evidently a 
consequence of the diffusive to ballistic motion conver- 
sion processes. Once a particle which comes from the 
LB is injected at x = x m id into a regular structure of 
the RB, it cannot become diffusive again and travels 
to x = £max ballistically. Moreover, it is unlikely for a 
diffusive particle to reach :r max , and indeed impossible 
if the length of the block tends to infinity, because the 
local current in the RB is negative. Hence almost every 
particle in the velocity regime v < 3.0 in Fig. [5] a) is a 
ballistic one. Note that the island structures can easily 
be identified with the ballistic islands in the RBs PSS 
(Fig. H]d)). In an analogous way we can understand the 
velocity regime v > 3.5: These particles are injected into 



FIG. 5. Properties of the escaping particles for the two block 
setup with the same parameters as in Fig(4] a) Phases and 
velocities at x = x max , c) at x = a: m j n . b) and d) show 
the corresponding normalised velocity distributions p(v). The 
inset in b) shows the velocity distribution for particles which 
arrive at x = £ m id for their first time. 



surface spanning curves above the RBs FISC at a; m id- 
Accordingly, they are not restricted to certain phases. 
The explanation why certain islands in the velocity 
regime 3.0 < v < 3.5 are avoided by particles at x max 
can be given straightforwardly after our previous discus- 
sions on the conversion process. As already mentioned 
in section IIV Al these islands correspond to regular 
structures of the PSS in the RB which have no overlap 
with the chaotic sea of the LBs phase space and hence 
no injection occurs. 

Finally, we turn our focus on the spreaded particles with 
3-0 < v < 3.5: A comparison with the corresponding 
PSS (Fig [U d)) reveals that these particles are located 
within the chaotic sea. Thus they have indeed passed 
the RB contrariwise to the directed current. To under- 
stand why this occurs predominantly in this velocity 
regime, a short detour on the typical length of Levy 
flights in the driven lattice is necessary. To this end we 
have simulated particles in a uniformly driven lattice 
with driving law d\(t) starting at x — with 2 ■ 10 6 
different initial conditions covering the phase space 
interval (0 < </>;„; < 27T, —5 < t>i n ; < 5). For every 
initial condition the number of barriers that the particle 
passes, before the sign of its velocity changes, is recorded 
and shown in Fig. [5] (initial conditions leading to 
regular motion were excluded and are shown in white). 
Apparently, the number of passed barriers before the 
velocity is reversed can differ by several orders of 
magnitude and strongly depends on the initial condition. 
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FIG. 6. Length of Levy flights for particles exhibiting chaotic 
dynamics in a uniformly driven lattice with driving law d\ (t) 
(parameters as in Fig. [2]d)). The white regions correspond 
to initial conditions leading to regular motion. 



Most interesting for our purpose is the observation that 
particles started in the velocity regime 3.0 < v < 3.5 
exhibit extraordinary long 'ballistic like' flights, which 
can be at the order of a few thousand barriers. Hence it 
is more likely for a particle -that remains diffusive once 
it passes x m id- to reach x max before being transported 
back to £ m id if it is within this region of extraordinary 
long Levy nights. Indeed, a comparison of Fig. |6] with 
Fig. [5] a) reveals, that the regions of long Levy flights 
coincide with the ones where diffusive particles reach 

Even though the overall appearance of the plot shown in 
Fig. [6] does strongly depend on the used parameters in 
the driving law, it is -for later usage- worth emphasising 
that the tendency for fast particles to exhibit much 
longer Levy flights than slower ones is a rather general 
feature in the driven lattice. This is mainly caused 
by two facts: Firstly, the lattice becomes a smaller 
perturbation for faster particles. Hence the average 
velocity change at a collision with a barrier is small for a 
particle which is close to the FISC (note that it is indeed 
zero for particles on the FISC) and as a consequence it 
takes many collisions with the barrier before a notable 
impact on the particle velocity occurs. Secondly, for 
particles with large initial velocities it is likely to become 
sticky to the FISC which -according to the discussions 
in [23]- leads to long ballistic flights. 
Besides leaving the system at x max , there is also the 
possibility for a particle to leave the system at X m i n . 
Again, phase and velocity at this particular position are 
recorded and shown in Fig. [5] c). This plot features 
qualitatively the same occupied domains as the one for 
i max but now mirrored at v = 0: Distinguished island 



structures belonging to a regular ballistic dynamics for 
less negative velocities, particles obeying chaotic dy- 
namics which avoid ballistic islands for —3.5 < v < —3.0 
and particles on regular spanning curves with velocities 
v < —3.5. The main difference appears to be the larger 
amount of spreaded diffusive particles for velocities with 
v ^ — 3 which do not correspond to ballistic islands. 
Both the difference as well as the similarity to the x max 
plot can be understood intuitively: Since the ensemble 
is initially located around the center of the LB, the 
number of barriers that these particles have to pass 
diffusively and opposite to the direction of the local 
current is roughly 5000, while it would be 10000 for 
particles that have reached X m id- Consequently, most 
of the spread diffusive particles seen in Fig. [5] c) are 
particles which have never reached x m id- In contrast to 
this, the particles reaching x m ; n within either ballistic 
islands or regular spanning curves have obviously passed 
^mid at least twice, because otherwise they could not be 
injected into the corresponding regular structures of the 
LB. Note that once a particle reaches a; m td, it is injected 
into the PSS of the RB with d\ (t) for positive velocities 
and into the PSS of the LB with do(t) for negative 
velocities. Since we chose do(t) = —di(t), the ballistic 
islands apparent in Fig. [5] c) are the same as in Fig. [5] 
a), but mirrored at v = 0. 



D. Injection probabilities into different regular 
structures 



So far we have seen that the two block system allows 
for diffusive to ballistic motion conversion processes. 
In this section we further investigate this phenomenon 
and discuss how likely injections into different regular 
regimes such as ballistic islands or spanning curves 
above the FISC are. 

To get some insight it is instructive to compare the 
normalised phase integrated velocity distributions at 
£max (Fig. [5]b)) with the one at £ m i n (Fig. E]d)). The 
for our purpose crucial observation is that the peak at 
high velocities (\v\ > 3) is less pronounced at x m [ n , and 
that the peaks corresponding to the ballistic islands 
appear slightly stronger populated (hardly visible) 
compared to the peaks in the distribution at x max . This 
effect can partially be explained by the larger number 
of diffusive particles at x min , leading to a broadening 
of the high velocity peak, but is also caused by certain 
characteristics of the injection process at x m jd as we 
shall explain in the following. 

A further understanding can be obtained by considering 
the velocity distributions for outgoing particles for a 
given number of times a particle has crossed x m id before 
it leaves the system (in the following we refer to this 
number as n CI ) as shown in Fig. [7] Apparently, for 
ri cr = 1 (Fig. [7] a)) the distribution features a very 
pronounced peak at velocities between 3 and 4, which 
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as in Fig. [U 



appears stronger populated compared to the peak in 
the total velocity distribution at x max shown in Fig. 
[5] b). On the contrary, the peaks corresponding to 
ballistic islands at smaller velocities are clearly weaker 
pronounced. In the n cr = 2 case it is less apparent, but 
the islands at small velocities are still less pronounced 
as for the total distribution in Fig. [5]d). For n cr = 3 
this effect is reversed and the low velocity peaks contain 
relatively more particles than they do in the total ve- 
locity distribution. Finally, particles leaving the system 
with n cr = 4 are very similar to the ones with n cr = 3 
and no obvious deviation in the velocity distributions is 
observed. It will become clear later, that this behaviour 
is strongly correlated with a different effect which is 
worth mentioning at this point: Not only the normalised 
distributions differ for different n cr , but the probabilities 
for a particle to leave the system after a certain number 
of crossing n cr are different as well. Hence, the injection 
probabilities into any kind of regular structures p Ucl has 
to be different for different n cr and we indeed found 
numerically, that the probability for a particle to be 
injected into a regular structure while it passes £ m id for 
the first time is p\ « 0.49, while it is P2 ~ 0.24 for the 
second passing of x m id and approximately pt ~ 0.15 for 
n cr > 2. 

To understand both the different overall injection 
probabilities p Uci for different n cl - as well as the different 
appearances of the corresponding velocity distributions, 
an argument which combines the length of Levy flights in 
different regions of the phase space (as discussed before 



and shown in Fig. [6] for the RB and to be mirrored at 
v = for the LB) together with the ergodicity property 
is required: Due to ergodicity parts of the phase space 
corresponding to long ballistic flights must be visited less 
frequently, but ones a particle gets there, it stays for a 
comparably long time. Hence, a particle initially started 
very close to x m id with a small velocity, reaches the 
interface after only a few collisions and as a consequence 
it is unlikely to reach the high velocity regime for such 
a particle. On the contrary, for a particle that started 
far -say some thousand barriers- away from x m id, it 
is likely that the particle reaches this high velocity 
domain at some point. Once it possesses such a high 
velocity, the length of its Levy flight is of the same 
order as its distance from ir m id and the particle typically 
reaches the interface while being still confined to this 
domain of phase space. We remark that this effect can 
easily be observed in the velocity distribution at ic m id 
for particles which reach this position for their first 
time as shown in the inset of Fig. [5]b). Apparently 
the distribution reveals a strongly pronounced peak at 
velocities close the LBs FISC in agreement with the 
previous discussion. Since these fast particles in the LB 
have velocities v > 3.5 they are injected into regular 
spanning curves above the FISC of the RB, which 
explains why particles at the first injection process have 
extraordinary high velocities and consequently why p\ is 
extraordinary large. Additionally, it illustrates why this 
domain appears to be strongly populated for n cr = 1 (cf. 
Fig. Ela)). 

Following the same arguments, the velocity distribution 
for n CI = 2 as well as the slightly enhanced value of 
p 2 compared to pi with i > 2 can be understood too. 
Because the particles are extraordinary fast at their 
first arrival at x m id, a comparably large fraction of 
particles which remain diffusive is injected into the 
part of the RBs PSS corresponding to long Levy flights 
(cf. Fig. [S]). Hence they surpass a large number of 
barriers in the RB which is at the order of 10 3 before 
their velocity is reversed for the first time. Afterwards, 
they are most likely transported towards x m id due to 
the negatively directed current in the RB. Now the 
same argument as before holds: because these particles 
are relatively far away from £ m id, they reach x m id 
predominantly with a high velocity. Hence, this domain 
is still strongly populated for the n cr = 2 injection 
process and pi is slightly enhanced. By now most 
of the fast particle have left the system and the still 
diffusive ones are transported back to x m id after only a 
few collisions. Consequently, the injection probability 
into the high velocity domain is suppressed for larger n CI . 



V. CONTROLLING THE VELOCITY 
DISTRIBUTION 

We have seen in section IIVI that a BL consisting of 
differently driven blocks allows for the conversion of 
diffusive- to ballistic motion and vice versa at inter- 
faces where the driving law changes. As a hallmark of 
these processes, we obtained velocity distributions with 
pronounced peaks at the velocities associated to regular 
structures in the underlying phase space. In the follow- 
ing section we investigate to what extend these conver- 
sion processes can be exploited to modulate the velocity 
distribution for outgoing particles in a controlled manner 
by adjusting parameters in the driving law. 
Again we focus on the two block setup (JVbi — 2, N — 
10 ) with two different driving laws do(t) for x < x m id 
and d\{t) for x > x m id- The parameters in do(t) re- 
main as before throughout the entire section and thus a 
positively directed current is induced in the LB. On the 
contrary, each parameter in di(t) is varied separately and 
its influence on the phase velocity distributions at x max 
and cc m in shall be studied. To this end we simulate an 
ensemble of 10 5 particles (initial conditions as in section 
IIV[) until all particles left the two block system. 
However, before we begin a detailed analysis of the phase 
velocity distributions for outgoing particles at x max in 
various parameter regimes it is sensible to have a brief 
discussion on their expected parameter dependence: As- 
suming a negatively directed current in the RB (which 
will be the case in most of the studied scenarios) , we can 
expect that the particles arrive at x max predominantly 
within either ballistic islands or regular curves above the 
FISC corresponding to the RBs phase space. Hence the 
phase velocity distributions at x max provide an 'image' 
of all the regular structures in the PSS of the RB which 
have overlap with the chaotic sea of the LBs PSS (Fig. 
0c)). Additionally to the ballistic particles some diffu- 
sive particles are expected to reach x max if the underly- 
ing phase space of the RB possesses regions of long Levy 
flights. Provided with these arguments we are able to ex- 
plain most of the phenomena occurring in the following 
sections. 



A. Frequency variations 

To begin with, we explore how a change of uj\ man- 
ifests itself in the particles' velocity distributions when 
they leave the system at either a; max or x m - m . The other 
parameters of d\(t) are kept constant at A\ w 0.57 and 
A$i =0 (inducing negatively directed currents in the 
RB for all values of cji). 

Figs. [5] a)-d) show the phase velocity distribution at 
i max with varying frequencies u>\ = 0.05, 1.50, 1.65 
and 4.00. Fig. [8] a) corresponds to the low frequency 
(cji = 0.05) regime where the barriers within the RB 
move very slowly compared to the barriers inside the LB. 
It shows a broad velocity band which is sharply confined 
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FIG. 8. Phase velocity distributions at x = £ max for a two 
block system with do(t) as in Fig. [4] The orange line indicates 
the FISC in the PSS of the RB. For di(t) the parameters are 
A x » 0.57, A$i = and ui = 0.05, 1.50, 1.65, 4.00 for 
a)-d). For e)-h) we have A$i = 0, wi = 1.0, and A\ = 
0.02, 0.40, 0.50, 1.00. For i)-l) we have Ai w 0.57, Wi = 1.0, 
and A$i = O.Itt, 0.2tt, 0.5tt, I.Ott. 



between v « 1.4 and v ss 3.8 where particles reach a; max 
without any apparent restrictions on their phases. This 
can be explained straightforwardly, by considering the 
expansion of the chaotic sea. As argued in [23( , the ve- 
locity for particles on the first invariant spanning curve 
(FISC) limits the chaotic sea and can be estimated by 
v± = ± d^ nax ± \J2Vq where ± correspond to positive 
or negative velocities and d^ lax is the maximal barrier 
velocity in either positive- or negative direction. For a 
small frequency u)\ the maximal velocity of the barrier in 
any given direction becomes small as well, leading to a 
FISC in the RB at very low velocities (for illustration we 
show the FISC in the RBs PSS as an orange line in FigE]). 
According to our previous discussion, almost all particles 
are injected into regular spanning curves above the RBs 
FISC at a; m id and are able to reach x max at arbitrary 
phases. The sharp cutoff at v ra 1.4 can be explained by 
employing the approximation of a static potential which 
yields a minimal particle velocity of v = \/2Vq « 1.41 to 



surpass the barrier. 

Increasing uj\ in the RB leads us to a second regime, 
where the FISCs in both blocks are at similar velocities 
and we show the corresponding phase velocity distribu- 
tion for cji = 1.50 in Fig. |8]b). The distribution reveals 
the expected domains of particles in ballistic islands and 
some diffusive particles close to the FISC in the RBs 
phase space. Both features are straightforward conse- 
quences of the previous discussions and we refrain from 
a reiteration of the arguments. However it is appealing 
to explore how the overall picture changes for a small 
variation in u>\. To this end the phase velocity plot for 
another frequency of u>x = 1.65 is presented in Fig. [8]c). 
Apparently, the phase velocity distribution of particles at 
i mM looks similar to the one for u)\ = 1.50. Nevertheless 
there are some differences worth emphasising. First of 
all, the 'large' islands survive the frequency change, but 
are shifted to slightly higher velocities for an increased 
lj\. Additionally, the islands at v « 4.0 appear now emp- 
tied around their center. By comparing the PSS of the 




FIG. 9. Phase-integrated velocity distributions at a; max (setup 
as as in Fig. [4| on a logarithmic scale for variations of a) uui , 
b) A\ and c) A$i. Remaining parameters as in Fig. [8] 

LB and the RB we see that those depleted parts of the 
ballistic islands in the RBs phase space are at velocities 
which are partially above the FISC of the LB. Since the 
injection at x m id can only occur for regular curves within 
the island that are at least in parts covered by the chaotic 
sea of the LBs phase space and the outer curves in each 
island reach to lower velocities than the inner curves, the 
latter ones are depopulated first. 



Finally, we address the regime where the FISC in the 
RB is at much higher velocities than the FISC in the 
LB, which is realised for uj\ >> luq and the phase ve- 
locity distribution is shown exemplarily for uj = 4.0 
in Fig [8] d). Apparently, only the island initially at 
(</) w 37r/2, v ~ 1.2) survives, but is shifted to consid- 
erably higher velocities (v w 3.0). Moreover, we notice 
a comparably large portion of diffusive particles at high 
velocities (7.0 < v < 9.0) below the FISC in the RB. 
To visualise the results over a broad range of uj\ , we in- 
tegrate the phase velocity distributions over the phase 
and normalise each in the range u)\ — 0.05 — 7.00 
(see Fig[9]a)). Evidently, the particles arrive within a 
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FIG. 10. Phase velocity distributions for particles which exit 
at :r m i n for a) uj\ — 0.05, b) uj\ = 1.50 and c) uj\ = 4.00 
(remaining parameters as in Fig(4]). The orange line indicates 
the FISC position in the LB. 

broad range of velocities for small wi's, while they are 
restricted to certain comparatively narrow velocity in- 
tervals at higher frequencies. Both can be understood 
within the above analysis: For small frequencies the par- 
ticles are predominantly injected into regular spanning 
curves and for higher ones mainly into individual ballistic 
islands. These islands are shifted to higher velocities for 
increasing ui and once they pass the FISC of the LB they 
simply disappear. What these velocity distributions re- 
veal additionally is that the islands mean velocities tend 
to increase linearly with wi. This is caused by the fact, 
that ballistic islands correspond to trajectories synchro- 
nised with the barrier oscillations in a way that every 
collision with the barrier occurs at distinguished phases. 
Each island is thereby characterised by its winding num- 
ber n which is defined (within a block) as the number of 
unit cells the particle passes within one period T. Hence 
the average velocity of a particle trapped in such an is- 
land is given by v n = 4? = 4^ • uii and thus proportional 
to uj\. 

Besides the possibility for leaving the system at x max 
particles can also exit at a; m i n . Although, the parame- 
ters in the LB are kept constant, the parameters in d\(t) 
can have a substantial influence on the particle dynamics 
at x m i n due to conversion processes after multiple cross- 
ing of x m jd . This is demonstrated in Fig. [TO] where the 
phase velocity distributions are shown exemplarily for 
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FIG. 11. Transmission probability pt as a function of a) wi, 
b) Ai and c) A^i. Remaining parameters as in Fig. [5] 



Wi = 0.05, 1.50 and 4.00. Again let us discuss first the 
small frequency regime (Fig. [TO] a)). Since the FISC in 
the RB is at very low velocities for u\ = 0.05, all the fast 
particles are immediately injected into regular curves and 
therefore leave the system at x = x max . Consequently, 
there is little probability for a particle to be re-injected 
into regular structures in the LB once it has entered the 
RB and as a consequence all the particles seen in Fig. 
[TO] a) are diffusive ones. More precisely, these are mainly 
particles that never reached x m id, i.e. the 'diffusive back- 
ground' which is always present at x m i n , independent of 
the parameters of di(t) and consists of approximately 
20% of all particles (note that this number tends to zero 
if the size of a block tends to infinity) . 
For increasing u>i the situation changes substantially, as 
demonstrated for u>\ — 1.50 in Fig. [TOJb). Evidently, 
some of the ballistic islands at minor negative velocities 
are filled. Additionally, we observe particles on regular 
spanning curves below the FISC, which were not present 
in the case U\ — 0.05. However, before we give an expla- 
nation by means of the underlying PSS we remark that 
the relevant parts of the phase space are now the ones for 
negative velocities (cf. Figs. [2]c) and d) for v < 0). By 
keeping this in mind, the difference of Fig. [TO] b) com- 
pared to Fig. [TO] a) can be understood intuitively by the 
FISC position: Since the FISC for negative velocities in 
the RB is now lower than the one for negative velocities 
in the LB, there is significant overlap of the chaotic sea in 
the RB and ballistic islands- as well as invariant curves 
below the FISC in the LB. Consequently, these parts of 
phase space in the LB can be populated by diffusive par- 
ticles in the chaotic sea of the RB which are transported 
back to a; m id by the negatively directed current. 
The same arguments hold for the high frequency domain 



that is exemplarily shown for u>\ = 4.00 in Fig. [TO1 c). 
Since the FISC for negative velocities in the RB is at 
lower, i.e. more negative, velocities compared to the case 
u>i = 1.50 the particles which penetrate into the RB dif- 
fusively can now be injected into regular curves of the 
LB at more negative velocities and therefore reach x m ;d 
faster than before. Note that even though the particles 
have no apparent modulation in phase, they are not uni- 
formly distributed within the accessible range of veloci- 
ties, which can be traced back to the underlying PSSs. 
For example the local minimum at v ~ —6.5 (cf. Fig. [TO] 
c)) is caused by a chain of ballistic islands embedded in 
the RBs chaotic sea which is decreasing the overlap and 
therefore the injection probability into regular structures 
of the LB within this velocity regime. 
We conclude this chapter with some remarks on the im- 
pact of a varying frequency in the RB on the transmission 
probability pr through the interface, i.e. the ratio of the 
particles reaching x max and the total number of particles 
reaching x m id- Its dependence on w\ is shown in Fig. 
[TT1 a) (note that we omit the diffusive background, i.e. 
we only take particles into account that have actually 
reached x m id at least once) . For small frequencies almost 
all particles are transmitted, while this changes drasti- 
cally for u>i > 1.0 when px starts to oscillate around 
Pt ~ 0.2 until it is again decreased above uii ~ 4.5 to 
Pt ~ 0.05 where it appears to saturate. All three regimes 
can be understood by means of our previous discussions: 
Since the particles arriving at x m ;d for their first time 
have predominantly high velocities (3.5 < v < 4.0) (cf. 
the inset of Fig. [S]c)), p\ (and hence px) is large when- 
ever regular structures in the phase space of the RB are 
located within this velocity regime. For small wi this 
is the case because the PSS in the RB is in this velocity 
regime filled with regular curves above the corresponding 
FISC. For wi > 1 this is no longer true (cf. for example 
[5] b)) and thus the particles rely on injection processes 
into ballistic islands rather than regular spanning curves 
to surpass the RB ballistically. Consequently pt drops 
significantly. For larger frequencies p\ (and therefore pt) 
increases whenever one of the ballistic islands in the RBs 
phase space is in the velocity regime 3.5 < v < 4.0 (cf. 
Fig. [9] a)) and decreases once this island disappears be- 
cause it is above the FISC of the LB, which explains 
the oscillatory behaviour of pt for 1.0 < Wi < 4.5. At 
uj\ ~ 4.5 the last ballistic island corresponding to an in- 
teger value of the winding number n disappears and thus 
Pt drops further. 



B. Amplitude variations 

Analogous to the previous discussion, we explore how 
a change of the oscillation amplitude A\ in the RB affects 
the particle properties for fixed values of Wi = 1.0 and 
A$i = 0. As before the parameters in do(t) remain un- 
altered. The corresponding phase velocity distributions 
at x max are shown exemplarily in Figs. |S]e)-h). The 
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regime of small amplitudes is represented by Fig. [8] e) 
which shows the particles phase velocity plot at x max for 
A\ = 0.02. Apparently, the particles are restricted to the 
same velocity interval as for the 'small frequency regime' 
illustrated in [8] a), namely 1.4 < v < 3.8. The interpreta- 
tion in terms of a quasi static barrier motion as provided 
previously for the case of small frequencies holds also 
here for the case of a small amplitude motion. Accord- 
ingly, most particles are injected into regular spanning 
curves above the FISC of the RB once they pass x m id 
for the first time. If we increase Ax the FISC is shifted 
to higher velocities and reaches v w 3.0 for A\ = 0.40. 
The associated phase velocity plot is shown in Fig. [8] 
f). For velocities below the FISC we observe the famil- 
iar island like structure. Particles can still be injected 
into curves above the FISC, where a chain of islands is 
avoided (v ~ 3.2). The latter one is caused by the fact, 
that these ballistic islands in the phase space of the RB 
have no overlap with the chaotic sea in the LB. This is 
due a chain of larger ballistic islands in the LBs phase 
space (Fig. c)) at the same velocity. 
Fig. [5]g) shows the corresponding graph for a slightly in- 
creased amplitude A\ = 0.50 and reveals that the FISC is 
-as expected- shifted to higher velocities. Furthermore, 
a region containing a considerable number of diffusive 
particles evolves at 3.0 < v < 3.5 followed by an island 
structure for even lower velocities. The surprisingly high 
number of diffusive particles at velocities slightly below 
the FISC is caused by extraordinary long Levy flights 
within this region due to a chain of cantori at v ~ 3. As 
a consequence of the small flux through this chain, the 
dwell time for particles in this part of the phase space 
is enhanced drastically 23]. According to our previous 
discussions in section IIVCI this enhances the probability 
for diffusive particles to reach x max in this part of the 
phase space. 

Fig. [8] h) shows the phase velocity distribution at x max 
for a large amplitude A\ = 1.0. Evidently, the FISC in 
the RB is in this case at higher velocities than it is in 
the LB. Thus, mainly particles within ballistic islands 
reach a; ma x- As we observed for a high frequency (cf. 
Fig. [S] c )), some of the islands are partially above the 
LBs FISC and their inner curves are therefore depleted. 
The phase integrated results over a range of amplitudes 
A\ = 0.02 — 1.0 is shown in Fig. Ob). Similar to what we 
observe for the frequency dependence the particles cover 
a broad velocity interval (1.5 < v < 4.0) for small ampli- 
tudes. With increasing A\ this interval decreases and ad- 
ditional narrow velocity intervals emerge. The latter are 
due to ballistic islands in the PSS of the RB. We remark 
that opposite to the frequency dependence (Fig. [9] a)), 
the corresponding velocity peaks are barely affected by 
changes of the amplitude. At least qualitatively, this be- 
haviour can be understood by considering trajectories of 
particles within these ballistic islands: On the one hand 
the constant mean velocity of an island can be under- 
stood by remembering that ballistic islands correspond 
to synchronised orbits, where the particle collides with 



the barrier at distinguished phases. Hence, their mean 
velocity only depends on the frequency rather than on 
the amplitude of the oscillation. On the other hand the 
precise point of appearance as well as the shape of the 
islands can very well depend on A\ and has to be deter- 
mined by numerical simulations. 

The amplitude dependence of the transmission proba- 
bility pt is shown in Fig. [11] b), where -as before- we 
consider only particles that reach a; m id at least once. Ev- 
idently, pt remains approximately unity for < A\ < 0.4 
and thereafter decreases steadily with further increas- 
ing A\, besides a weakly pronounced local maximum at 
A s=s 0.8. Following the arguments presented in the pre- 
vious section, this behaviour is a result of the overlap of 
the chaotic sea in the LB and regular structures within 
the RB. 

No relevant new phenomena are observed for the particles 
exiting at ic m i n which is why we refrain from providing a 
discussion of this case. 



C. Phase variations 

Let us finally explore the impact of phase changes A$i 
for fixed u>i = 1.0 and A\ m 0.57. The phase velocity dis- 
tributions at £ max are shown in Fig. |5]i)-l). Fig. [5]i) and 
j) correspond to comparably small phase differences of 
A$i = 0.1-7T and A$i = 0.2-7T respectively. A compar- 
ison of these graphs with the one for A$i = (Fig. [5] 
a)) reveals that the main difference is a small shift of 
the entire phase space to the region of smaller phase val- 
ues. Consequently, the n — 4 island (v sa 3.2) in the 
RB is not completely covered by the corresponding is- 
land in the LB anymore. Hence, the outer most curves 
can be populated by particles. The tendency that islands 
are moved to smaller values of the phases for increasing 
A$i is still apparent for A$i = 0.5-7T as seen in Fig. 
[5] k) . Additionally, the n = 4 island is fully populated 
by particles and more diffusive particles manage to reach 
i max . Finally, there are no particles in spanning curves 
above the FISC anymore. In contrast to the frequency 
and amplitude dependence, A$i has a substantial in- 
fluence on the symmetries of the Hamiltonian and thus 
on the transport (cf. FigJ3]). In fact for A$i = 0.57T 
time reversal symmetry is restored (dx(t, A<I>i = 0.57r) = 
rfi(— t, A$i = 0.57!")) and the directed current vanishes. 
Obviously, this increases the probability for a particle to 
traverse diffusively the RB, which explains the notable 
amount of non ballistic particles in Fig. [5]k). Finally, 
the transport in the RB is reversed, i.e. points in a pos- 
itive direction, for A$i = ix. Moreover we notice that 
di(t, A$ = 7r) = — di(t + 7r, A$ = 0), i.e. the driving law 
in the RB equals the driving law in the left one besides 
an initial phase shift of tt. Consequently, the ballistic is- 
lands for d\(t) are at the same velocities but at phases 
shifted by tt. Since the PSS for doit) is 'almost' invariant 
under a shift of tt, the overlap from the chaotic sea in the 
LB and regular structures in the right one is comparably 
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small for A$i = tt. Thus, very few particles reach x max 
ballistically and the phase velocity distribution is domi- 
nated by chaotic particles (Fig. [8]1)). 
Fig. |H]c) shows the A$i dependence of the phase inte- 
grated velocity distributions. We observe an increasing 
amount of diffusive particles for A$i close to A$i = rr 
compared to A$i = (or A$i = 2w). As stated before, 
this is a consequence of the dependence of the direction 
of the current in the RB which is positively directed for 
A$! = 7r while it is negatively directed for A$! =0. 
As a last remark on the two block setup the transmission 
probability pr as a function of A$i is shown in Fig. \TT\c) 
and reveals an oscillatory behaviour with local maxima at 
A$i = 0, 7r and 2i\. For frequency and amplitude varia- 
tions we argued that the value of pt is determined by the 
overlap of regular structures in the RB with the chaotic 
sea of the LB. However, in the present case, this over- 
lap is minimal for A$i = n. Apparently, this decreasing 
overlap for A$i — » tt is compensated by an increasing 
probability for a particle to surpass the RB diffusively 
due to the positively directed current (cf. FigJ3]). 



VI. VELOCITY DISTRIBUTIONS IN 
SUPERLATTICES 



In the previous sections we have demonstrated how 
a setup build up out of two blocks with different driv- 
ing laws allows for conversion processes from diffusive- 
to ballistic motion. Even more we were able to control 
the velocity distributions for outgoing particles at x max 
by adjusting parameters in the RB. In the following we 
argue how the so far gained insights can be exploited 
to maintain mono energetic- pulsed particle beams out 
of diffusive particle ensembles in superlattices containing 
a few hundred blocks. The general outline of the used 
scheme is as follows: We start with an initially diffusive 
particle ensemble in the B = block (with B being the 
block index, cf. Fig [1]) which is transported towards a 
first interface where particles can be injected into ballis- 
tic islands of the B — \ blocks phase space. The param- 
eters in the driving laws are chosen such that these now 
ballistic particles travel opposite to the directed currents 
and thus we obtain a peaked velocity distribution at the 
end of the B = 1 block. For the following blocks with 
B = 2, ..., 100 we show how the width of each peak in the 
velocity distribution can be tuned by adjusting the am- 
plitude of the barrier oscillation blockewisely. As a last 
step, we demonstrate how an appropriate choice in the 
driving laws for B > 100 allows to preserve one of the 
peaks in the velocity distributions while the other peaks 
are subsequently removed. Thus we obtain a mono en- 
ergetic particle beam for outgoing particles in the super- 
lattice. Moreover, the beam is pulsed in a sense that the 
particles leave the systems only at distinguished phases. 



A. Interface dynamics of ballistic particles 

Before we start a detailed discussion of the physics in 
the BL containing a few hundred blocks, let us again con- 
sider the simple case of a two block setup to introduce a 
new type of conversion processes which occurs in larger 
BL namely ballistic to ballistic- or ballistic to diffusive 
conversion. To make our discussion more concrete we 
consider again a setup with driving laws as in section HVl 
Hence the PSS for the LB is shown in Fig. [5] c) and 
the PSS for the RB is shown in Fig. [2]d). In contrast 
to the previous discussions we explore the possible con- 
version processes for a ballistic particle arriving at the 
interface. For example, consider a particle beam started 
at %mm which uniformly occupies the ballistic island at 
(v = 1.8, (f> = 3n/2) (in Fig. H]c)) and passes ballistically 
the LB. Apparently, these particles would be entirely in- 
jected into the chaotic sea of the PSS in the RB once they 
pass a; m id. Due to the negatively directed transport in the 
RB the particles are transported back to a; m id where they 
can again be injected into any regular structures of the 
PSS corresponding to the LB, or after several passings of 
a^mid into regular structures of the RB. Hence, the out- 
going particles at x max for this initially mono energetic 
beam would occupy all accessible regular structures of 
the PSS of the RB and the corresponding velocity distri- 
bution at a; max would contain multiple peaks. 
As a second example we consider a particle beam (again 
started at £ m m) which passes the LB by uniformly occu- 
pying the ballistic island at (v — 3.2,0 = 3ir/2). In this 
case, some particles are injected into the ballistic island in 
the RBs PSS at similar coordinates, while others become 
diffusive. The particles which remain ballistic traverse 
the RB and cause a dominant peak in the velocity distri- 
bution at x max . For the particles which become diffusive 
the same arguments hold as before. Thus these particles 
lead to less pronounced peaks in the velocity distribution 
at velocities corresponding to any kind of regular struc- 
ture in the RB. Accordingly the initial particle beam was 
converted into a particle beam with a smaller width, be- 
cause the ballistic island in the LB, i.e. for the initial 
beam, is larger than the island in the RB in which these 
particles arc injected. In addition to this, peaks in the 
velocity distribution emerge due to injection of diffusive 
particles after multiple crossings of x m id- 



B. Amplitude variations in superlattices 

In the following section we demonstrate how the pre- 
viously discussed interface dynamics of ballistic particles 
can be exploited to narrow the velocity distribution of 
particles in an appropriately designed superlattice. To 
this end we consider a setup build up out of Nbi = 101 
blocks (whereas each block contains N — 10 4 barri- 

= to x max = NLNsi- 



ers) which expands from x n 

Accordingly, the positions of the interfaces, i.e. the 

positions where the driving laws change, are given by 
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x B = NLB with B = 1, ..., 100. For the B = block 
the driving law is d (t) with parameters as in Fig. [2]c) 
inducing a positively directed current. For B = 1, .., 100 
the driving laws are: d B (t) = A B [cos{2.2t) + sin(4.4t)] 
with A B = 0.3 + 0.07 • B (i.e. the amplitude is slowly 
increased from 0.3 to 1.0) inducing negatively directed 
currents. The initial conditions for the simulated ensem- 
ble are t = 0, 0ANL < x < 0.6NL and -0.1 < v < 0.1. 
Hence the particles are located within the chaotic sea of 




FIG. 12. Velocity distributions at positions a) x% = 2NL and 
b) a; max = 101 NL. Parameters in do(t) as in Fig. [4] For < 
B < 100 we set d B (t) — A B [cos(ujB't + 'PB)+sm(2(ijjB't + ^>B))] 
withes = 2.2, <p B =0 and A s = 0.3 + 0.07(5-1). The inset 
in a) shows p(y) at X = 200NL with ljb — 2.2, Ab = Aioi 
and ip B = tv{B - 1) for 100 < B < 200. The inset in b) 
shows p(v) at x — 500NL with ujb — 2.2, Ab = -Aioi and 
ip B = 0.02tt(B - 1) for 100 < B < 500. 



peak (v « 1.8) is associated to a n = 1 island. 
In the following blocks the amplitude of the barrier os- 
cillation is subsequently increased and as a result we ob- 
tain a velocity distribution with two narrow peaks at the 
same velocities as before for particles at a; max = 101 NL 
(Fig. HUb)). Additionally, we observe some particles 
with velocities v > 4.5. The reason for the two dominant 
peaks is that the amplitude has (as argued in section |V| 
only little influence on the position of ballistic islands 
in phase space. Thus, most particles remain ballistic at 
each interface. However the amplitude does have a no- 
table influence on the size of the islands and by choos- 
ing the amplitude appropriately, one can tune the width 
of the velocity distribution by adjusting the size of the 
corresponding ballistic islands. In the present setup we 
exploit that an increasing amplitude leads to a decreas- 
ing size of the islands for the used parameters. Hence, 
the velocity distribution is squeezed when the particles 
propagate further into the super lattice. The fast particles 
with v > 4.5 correspond to particles in ballistic islands 
of the underlying phase space which is best seen in the 
phase velocity distribution at x = 101NL (Fig. [T3"lb)b 
In fact, these are the in the previous section discussed 
peaks in the velocity distribution that emerge due to in- 
jection of diffusive particles after multiple crossings of an 
interface. 



C. Peaked velocity distributions in superlattices 



the B — block with driving law do(t) and transported 
towards the first interface at xi — NL. At this point they 
can be injected into ballistic islands into the phase space 
in the B = 1 block. Since the local current in this block 
is negatively directed, it is hard to surpass for diffusive 
particles and we obtain a peaked velocity distribution at 




FIG. 13. Phase velocity distributions at a) x — 2NL and b) 
X = 1Q1NL. Orange line indicates the FISC. (Parameters as 
in Fig. El). 

x 2 = 2NL (Fig. [T2]a)) which is is dominated by a peak 
at v « 3.5 and a less pronounced one at v « 1.8. The 
phase velocity distribution at x 2 — 2NL (Fig. QSla)) re- 
veals that the dominant peak (v sa 3.5) can be related to 
an island with winding number n = 2, while the second 



The last step to a mono energetic particle beam is to 
remove one of the peaks in Fig. [12] b) without losing too 
many particles in the other one. This can be done by ex- 
ploiting the symmetries of both islands, which is achieved 
by adding more blocks to the superlattice with driving 
laws: d B (t) = 1.0[cos(2.2i + <p B ) + sin(4.4i + 2.2tp B )} for 
101 < B < 201 with (p B = tt(B - 1). Before we show the 
resulting velocity distributions, let us briefly discuss the 
idea behind the chosen driving laws: On the one hand 
we have seen that the peak at v ~ 3.5 corresponds to 
a n = 2 island and consists of two island structures at 
the same velocity but at different phases (cf. Fig. IT3"1) . 
Hence, an additional phase shift ip = ir in the driving law 
'maps' both island into each other and most particles re- 
main ballistic. On the other hand the n = 1 island which 
is responsible for the peak at v fw 1.8 is 'mapped' into 
the chaotic sea for such a phase shift and particles in it 
become diffusive. Even though some of these now diffu- 
sive particles might be reinjected into a ballistic island 
of the following block, the majority is transported away. 
Consequently, after performing this procedure multiple 
times, one obtains a mono energetic particle beam. The 
resulting velocity distribution is shown in the inset of Fig. 
[T2l a) and reveals that we obtain indeed the desired form 
of a mono energetic particle beam. 

At this point we remark that the described technique 
of removing peaks according to the symmetry of their 
associated ballistic island works for a wide range of dif- 
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ferent parameter values as well as for ballistic islands 
with higher winding numbers. Unfortunately, it does 
not apply for the n = 1 island and thus we can not 
use it to remove the peak at v sw 3.5 while keeping 
the one at v s=s 1.8. However, we can exploit that for 
a large amplitude the n = 1 island tends to cover a 
larger range of phases in the PSS (cf. Figs. [5] g) and 
h)). Hence, a small phase shift in the driving law re- 
moves relatively fewer particles in the n = 1 island com- 
pared to the ones with higher n. Following this idea, we 
choose d B (t) = A 101 [cos(2.2t + tp B ) + sin(4.4t + 2.2ip B )} 
for 101 < B < 501 with ip B = 0.02n(B - 100) and the 
resulting velocity distribution at x = 500NL is shown in 
the inset of Fig. [T2] b). Again, we obtain the desired 
distribution of a monoenergetic particle beam. 



VII. CONCLUSION 

We have explored the classical non-equilibrium dynam- 
ics of particles in a one-dimensional driven superlattice 
which consists of blocks each containing many individual 
barriers. While similar systems that are usually studied 
in this context consist of lattices where all barriers are 
governed by the same time-dependent force, i.e. driv- 
ing law, we allowed for a different driving in each block. 
In doing so we show that the thus obtained variability 
leads to remarkable new dynamical phenomena. To this 
end we analysed in detail how the blockwise variation of 
the driving law gives rise to conversion processes from 
diffusive- to ballistic motion and vice versa at the in- 
terfaces, i.e. the positions in the superlattice where the 
driving law changes. The combination of directed trans- 
port and these conversion processes enabled us to obtain 



peaked velocity distributions in a simple system contain- 
ing only two blocks with different driving laws providing 
oppositely directed currents. Additionally, we observed 
strong correlations between the phases and velocities for 
the escaping particles even though the initial particle en- 
semble is of exclusively diffusive character. Even more, 
we found that the velocity distributions as well as the 
correlations can be modified in a controlled manner by 
adjusting parameters such as frequency or amplitude in 
the driving. Finally, we present a scheme for superlat- 
tices containing a few hundred blocks by witch a dif- 
fusive particle ensemble can be converted into a pulsed 
particle beam, whose mean energy and width in momen- 
tum space can be adjusted. Since this scheme mostly 
depends on simple symmetry arguments it is viable over 
a wide range of parameters. Thus it should be applica- 
ble to experimental setups, such as layered semiconduc- 
tor heterostructures with different AC drivings or even 
to cold atom experiments in which counter propagating 
laser beams can create a one-dimensional lattice poten- 
tial. By passing the laser beams through two acousto- 
optical modulators the desired AC drivings can be ob- 
tained. 

As a future perspective it would be intriguing to explore 
-both theoretically as well as experimentally- the ana- 
logues of the presented effects in the quantum regime. 
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